Analyzing  the  Stability  of 
Floating  Ice  Floes 

Barry  Couternnarsh  and  Randy  McGilvary  December  1 994 


Abstract 

This  report  describes  an  experimental  study  to  measure  the  pressure  caused 
byfiuid  acceieration  beneath  a  fioating  parallelepiped  biock.  Dynamic  fluid 
pressure  was  measured  at  discrete  points  beneath  the  block  for  several 
flow  velocities,  flow  depths,  block  angles  of  atfack  and  block-thickness-to- 
depth  ratios.  The  measured  pressures  were  used  to  calculate  block  underturning 
moments,  and  a  hydrostatic  analysis  was  used  to  calculate  a  block  righting 
moment.  From  this,  a  densimetric  Froude  underturning  criterion  is  presented. 
The  measured  hydrodynamic  pressure  distribution  on  the  bottom  of  a  single 
model  ice  floe  is  used  to  estimate  the  dynamic  stability  at  three  thickness- 
to-depth  ratios.  The  energy-based  analysis  details  the  conditions  required 
for  instability,  metastability  and  stability.  At  three  thickness-to-depth  ratios, 
block  rotational  inertia  has  the  effect  of  reducing  the  Froude  stability  number 
by  5  to  1 0%  over  a  completely  static  stability  criterion. 
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Coven  Underturning  block  (I  =  block  length;  V  =  fluid  velocity;  9  =  rotation 
angle;  t  =  block  thickness;  H  =  flow  depth). 


For  conversion  of  SI  metric  units  to  U.S./British  customary  units  of  measurement 
consult  ASTM  Standard  E380-89a,  Standard  Practice  for  Use  of  the  International 
System  of  Units,  published  by  the  American  Society  for  Testing  and  Materials, 
1916  Race  St.,  Philadelphia,  Pa.  19103. 
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NOMENCLATURE 


A  block  plane  area 

C  dynamic  fluid  pressure  coefficient 

C(0*)  dynamic  fluid  pressure  coefficient  at  the  angle  of  maximum  hydrostatic  righting 
moment 

COPx  center  of  pressure  in  the  X-direction  (across  block  width) 

COPy  center  of  pressure  in  the  Y-direction  (along  block  length) 
f  force 

f j.  densimetiic  Froude  number 

Ff  critical  densimetric  Froude  number  for  block  stability 

g  acceleration  due  to  gravity 
H  flow  depth  upstream  from  the  block 
I  moment  of  inertia 
L  block  length 

M  moment  acting  on  the  block 
Mf  undertuming  moment  created  by  fluid  flow 

Mf(0)  imdertuming  moment  created  by  fluid  flow  at  an  angle  of  attack 

Mf(0*)  rmdertuming  moment  created  by  fluid  flow  at  angle  of  maximum  righting  moment 

Mj.  hydrostatic  righting  moment 
Mj.(0)  hydrostatic  righting  moment  at  an  angle  of  attack 
Mj.(0*)  hydrostatic  righting  moment  at  angle  of  maximum  righting  moment 
My  undertuming  moment  in  the  Y-direction 

Mx  undertuming  moment  about  X-axis 

My  undertuming  moment  about  Y-axis 

■^max  maximum  hydrostatic  righting  moment 
P(X,Y)  differential  pressure  from  hydrostatic  at  a  point  (X,Y) 

P  mean  differential  pressure  from  hydrostatic 

St  Strouhal  number 

t'  block  freeboard 

t  block  thickness 

V  fluid  velocity  upstream  from  the  block 
V*  approach  velocity  of  instability 

X  spatial  coordinate,  distance  across  block  width 

Y  spatial  coordinate,  distance  along  block  length 

Z  spatial  coordinate,  distance  perpendicular  to  X-Y  plane 

a  nondimensional  moment  coefficient  in  Y-direction 
a(0*)  nondimensional  moment  coefficient  in  Y-direction  at  maximum  righting  moment 
angle 

p  fluid  density 

p'  block  density 

pj  ice  density 

p^v  water  density 

0  block  angle  of  attack  (rotation)  from  horizontal 

0*  block  angle  of  attack  at  the  angle  of  maximum  righting  moment 

00  initial  block  angle 

0  block  angle  velocity 

0  block  angular  acceleration 
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INTRODUCTION 

Ice  floes  that  come  to  rest  against  an  ob¬ 
struction,  such  as  an  ice  boom  or  an  intact  ice 
sheet,  either  wfll  remain  there  or  undertum. 

Floes  that  do  not  undertum  become  part  of 
an  ice  cover  that  progresses  upstream,  while 
floes  that  undertum  may  become  lodged  be¬ 
neath  the  downstream  cover  and  contribute 
to  the  growth  of  an  ice  jam.  It  is  for  these  rea¬ 
sons  that  block  imdertuming  has  been  stud¬ 
ied  to  help  us  rmderstand  the  conditions  that 
affect  ice  floe  stability.  Figure  1  is  a  schematic 
diagram  of  the  physical  situation  that  defines 
important  variables. 

Of  the  studies  that  examined  block  stabil¬ 
ity  criteria,  most  have  assumed  simple  em¬ 
pirical  descriptions  for  the  forces  and  moments 
caused  by  fluid  acceleration  around  the  floating 
block.  In  this  study,  we  measured  the  actual  pres¬ 
sures  acting  on  the  bottom  surface  of  a  floating 
block,  considering  the  effects  of  flow  velocity,  flow 
depth  and  block  angle  of  attack  on  the  pressure  dis¬ 
tribution.  These  measured  pressures  were  both  posi¬ 
tive  (stabilizing)  and  negative  (destabilizing),  de¬ 
pending  primarily  upon  the  angle  of  attack.  The  re¬ 
sulting  integrated  force  magnitude  and  center  of 
pressure  locations  were  used  to  calculate  the  total 
resultant  imdertuming  moment  acting  on  the  block. 
To  simplify  the  experiment,  we  used  a  rigidly  sup¬ 
ported  block  to  eliminate  the  higher  order  dynamics 
associated  with  block  rotation,  such  as  the  fluid 
added  mass,  block  rotational  inertia  and  the  chang¬ 
ing  relative  velocity  between  the  fluid  and  the  block. 

We  compare  the  block's  buoyancy-induced 
righting  moment  as  a  function  of  attack  angle  with 
the  imdertuming  moment  to  determine  the  condi¬ 
tions  when  the  imdertuming  moment  would  ex¬ 
ceed  the  righting  moment  and  the  block  would 
undertum.  From  this,  we  propose  a  simple  stability 
criterion  that  relates  the  block  densimetric  Froude 
number  to  a  dimensionless  moment  coefficient  at 
the  angle  of  maximum  righting  m  oment. 


Figure  1.  Undertuming  block  problem. 

We  then  use  the  measured  hydrodynamic  pres¬ 
sures  to  estimate  the  dynamic  stability  at  three  thick- 
ness-to-depth  ratios.  This  energy-based  analysis  de¬ 
tails  the  conditions  required  for  instability,  metasta¬ 
bility  and  stability.  We  show  that,  at  all  three  thick- 
ness-to-depth  ratios,  block  rotational  inertia  reduces 
the  stability  Froude  number  by  5  to  10%  over  a  com¬ 
pletely  static  stability  criterion. 

BACKGROUND 
Literature  cited 

Past  investigators  have  assumed  a  uniformly  dis¬ 
tributed  undertuming  pressure  across  the  block  and 
conducted  static  analyses  on  data  that  included 
block  angular  acceleration  and  angular  velocity.  Per¬ 
haps  the  most  popular  stability  criterion  was  that  of 
Pariset  and  Hausser  (1961),  who  assumed  that  ice 
floes  would  become  unstable  when  the  upstream 
edge  lost  its  freeboard.  This  incipient  instability  con¬ 
dition  became  known  as  the  "no-spill  criterion." 

Uzuner  and  Kennedy  (1972)  developed  an  un¬ 
dertuming  criterion  where  block  rotation  was  the 
primary  mechanism.  Their  analysis  included  an  em¬ 
pirical  moment  coefficient  that  depended  upon  t/L 


and  p'/Pw  was  determined  by  disregarding  the 
effect  of  t/H  (where  f  is  the  block  thickness,  L  the 
block  length,  H  the  upstream  flow  depth,  p'  the 
block  density  and  p^  ^e  fluid  density).  They  also 
defined  a  "no-spfll"  condition  as  the  point  at  which 
the  upstream  top  comer  of  the  block  was  lowered  to 
the  stagnation  height  of  the  water  surface,  and  they 
declared  this  to  be  the  point  of  instability. 

Ashton  (1974)  employed  the  "no-spill"  criterion, 
finding  a  vital  dependence  upon  block  thickness 
and  block  thickness-to-depth  ratio.  But  he  believed 
that  thickness-to-length  ratio  was  of  little  impor¬ 
tance,  as  was  depth  alone,  although  he  noted  a  slight 
increase  in  instability  with  increasing  thickness-to- 
length  ratios.  He  also  provided  a  simplified  moment 
an^ysis,  taking  into  accoimt  the  block  weight,  glo¬ 
bal  fluid  acceleration  to  maintain  continuity  and  the 
displaced  volume  of  fluid.  In  his  analysis  the  hydro¬ 
static  righting  moment  was  proportional  to  the  sine 
of  the  angle  of  block  rotation,  with  the  moment  as¬ 
sumed  to  act  at  L/2  and  with  an  implicitly  assumed 
lift  coefficient  of  1. 

Coutermarsh  (1986)  studied  ice  floes  striking 
against  the  45°  sloped  face  of  a  military  floating 
bridge  to  see  if  they  would  jam  or  undertum.  That 
situation  was  very  similar  to  the  one  we're  con¬ 
cerned  with  here,  except  that  the  block  rotated  about 
its  upper  downstream  comer,  a  factor  that  made  a 
minimal  difference  in  the  analysis.  Tbiat  study,  how¬ 
ever,  pointed  out  that  the  undertuming  phenomenon 
was  poorly  described  for  the  deep  water  region  and 
that  no  one  previous  study  had  described  instability 
over  both  the  deep  and  shallow  water  regions. 

Daly  and  Axelson  (1990)  performed  detailed  cal¬ 
culations  to  determine  the  hydrostatic  righting  mo¬ 
ment  related  to  the  displaced  volume  of  the  block  as 
it  rotates.  They  found  this  righting  moment  to  be 
highest  at  an  angle  that  depended  on  block  density, 
thickness  and  length.  They  reasoned  that  the  block 
became  unstable  when  the  undertuming  moment 
caused  by  the  fluid  flow  was  greater  than  the  maxi¬ 
mum  righting  moment;  they  discounted  the  previ¬ 
ous  notion  of  the  "no-spill"  condition.  Their  analysis 
also  determined  a  moment  coefficient  that  included 
two  empirically  determined  parameters  found 
through  analysis  of  existing  laboratory  data. 


about  the  Y-axis  (transverse  to  the  stream  flow)  that 
would  be  calculated  as 

Mx  =  \AP{X,y)ydA  (1) 

MY  =  JAfl(XY)XdA  (2) 

where  P(X,Y)  is  the  pressure  differential  from  hy¬ 
drostatic  and  A  is  the  block  plane  area. 

With  the  total  force  at  the  bottom  surface  as 

F  =  JAP(X,Y)riA  (3) 

the  average  pressure  is 


II 

(4) 

We  define  a  dynamic  fluid  pressure 
(Euler  number  or  lift  coefficient)  as 

coefficient 

fpw 

(5) 

with  being  water  density  and  V  flow  velocity 

upstream  of  the  block.  The  centers  of  pressure  in  the 
X(COPx)  and  Y(COPy)  directions  are 

My 

COPv  =  — 

X  p 

(6) 

COPy=^. 

(7) 

In  our  study  we  were  only  concerned  with  the  Y- 
direction  moments  and  forces,  as  they  are  the  ones 
that  win  contribute  to  undertuming.  If  symmetry  is 
present  about  the  Y-axis,  the  X  moments  wfll  bal¬ 
ance,  i.e.,  the  block  wfll  not  tend  to  rotate  aroimd  the 
Y-axis.  We,  therefore,  define  a  nondimensional  mo¬ 
ment  coefficient  a  in  the  Y-direction  as 


a  =  C 


or  alternatively 


a  =  - 


-Pw  V^AL 


(8) 

(9) 


Analytical  discussion 

We  approach  our  block  instability  study  by  fol¬ 
lowing  past  convention  and  assume  that  the  block 
rotates  about  its  lower  downstream  comer  because 
of  a  pressure  reduction  at  the  bottom  surface  of  the 
block.  This  pressure  reduction  creates  moments 
about  the  X-axis  (parallel  to  the  stream  flow)  and 


where  L  =  floe  length 

C  =  lift  coefficient 

COP/L  =  center  of  pressure  normalized  against 
the  block  length. 

The  subscript  Y  has  been  dropped  from  the  COP 
here  (eq  8).  These  parameters  are  related  to  the  aver- 
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age  pressure  acting  on  the  block  underside  P  by  eq 
5  and 


copy  Mu 

L  J~ PAL 


(10) 


where  My  is  the  undertuming  moment  in  the  V-di- 
rection. 

In  addition  to  the  forces  working  to  undertum 
the  block,  there  is  a  hydrostatic  righting  moment 
that  depends  on  fluid  density,  block  density,  length, 
thickness  and  angle  of  rotation.  Daly  and  Axelson 
(1990)  showed  that  the  righting  moment  for  a  paral¬ 
lelepiped  block  is  represented  by 

^  - =  f '  sin(0)  (p'  - 1)  +  p'(sec  (0)  - 1) 

2Pw  giA 


(12c) 


02  =  cos-1  [p']  (12d) 

and  t  =  block  thickness 
Pi  =  ice  density 
g  =  acceleration  of  gravity 
Oj  =  angle  at  which  the  leading  edge  loses  its 
freeboard 

02  =  angle  at  which  the  trailing  edge  loses  its 
freeboard. 


[f '  tan  (0)  +  (f '  -  L  tan  (0))  sin  (0)  +  L  sec  (0)] 

+ y  tan^  (0)  -I-  CSC  (0)  -  y  sin  (0)  j  (11a) 

for  0  <  0  <  01  and 

^  — = (L  cos  (0)  -t  f  sin  (0))  (l  -  p') 

2pw  gt^ 

-  y  (l  -  p'  sec  (0))^  j^f '  +  y  (f  cos  (0) 

-P)(csc2(0)  +  l)j  (lib) 

for  01  <  0<  02  and 

Mr=y(pw-Pi)gf^ 


Daly  and  Axelson  calculated  this  moment  for  both 
floating  and  submerged  blocks,  but  the  form  of  the 
equations  they  presented  made  it  difficult  to  deter¬ 
mine  the  magnitude  and  reaction  point  of  the  right¬ 
ing  force.  To  determine  these  components,  we  cal¬ 
culated  the  hydrostatic  righting  moment  through 
the  use  of  hydrostatic  pressure  prisms  (Streeter  and 
Wylie  1979).  These  lead  to  the  nondimensional  hy¬ 
drostatic  moment  given  as  a  function  of  length  to 
thickness,  as  shown  in  Figure  2.  Note  that  for 
lengths  in  excess  of  about  five  thicknesses,  the 
maximum  hydrostatic  righting  moment  may  be 
approximated  by 

■^inax  ~  (1  ~  Pi/ Pw)  Pw  g^^  (b/Z) .  (13) 

In  addition,  the  angle  of  maximum  hydrostatic 
righting  moment  0*  is  shown  in  Figure  2,  and  for 
lengths  greater  than  about  2.5  thicknesses  may  be 
approximated  as 

0*  =  tan"^  [t/L] .  (14) 


(Lcos(0)-tfsin(0)) 

(11c) 

for02^0^7t/2 

where 

1! 

(12a) 

Pw 

(12b) 

We  have  also  plotted  the  previously  used  "no¬ 
spill"  angle  criterion  in  Figure  2  for  comparison.  It  is 
important  to  note  that  the  no-spiU  angle  is  signifi¬ 
cantly  smaller  than  the  angle  at  which  the  maxi¬ 
mum  righting  moment  occurs.  For  our  particular 
block,  the  hydrostatic  righting  moment  is  shown  in 
Figure  3  for  f  =  0.076  m  and  L  =  0.62  m.  Our  block- 
length-to-thickness  ratio  is  8  and  the  approxima¬ 
tions  to  the  maximum  righting  moment  and  angle 
of  maximum  righting  moment  are  in  very  good 
agreement  with  the  actual  values. 
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Figure  2.  Hydrostatic  righting  moment.  Both  the  calculated  and  approxi¬ 
mate  angle  at  which  the  maximum  hydrostatic  righting  moment  occurs  are 
plotted  vs.  the  block-length-to-thickness  ratio.  Our  block  L/t  =  8.  Also  shown 
is  the  previously  used  "no-spill"  angle  condition  for  block  instability. 


Figure  3.  Calculated  block  hydrostatic 
righting  moment  for  our  particular 
block  vs.  angle  of  attack.  Both  the  calcu¬ 
lated  and  approximate  angles  of  maximum 
righting  moment  are  also  shown.  The  hori¬ 
zontal  line  is  the  approximated  maximum 
righting  moment  for  our  block. 

EXPERIMENTAL  METHODS 
Objective 

Our  objective  was  to  measure  the  pressures  act¬ 
ing  on  the  block.  As  Daly  and  Axelson  (1990) 
pointed  out,  most  of  the  previous  work  has  implic¬ 


itly  assumed  that  the  flow  beneath  a  block  is  one¬ 
dimensional,  resulting  in  a  uniform  dynamic  pres¬ 
sure.  (Larsen  [1975]  assumed  a  two-dimensional,  in- 
viscid,  irrotational  flow  to  calculate  a  one-dimen¬ 
sional  pressure  variation  on  the  block  tmderside.) 
We  felt  that  three-dimensional  flow,  producing  the 
two-dimensional  pressure  distribution  on  the  block, 
better  represented  actual  conditions.  Furthermore, 
we  assumed  that  the  flow  pattern  would  be  prima¬ 
rily  influenced  by  flow  depth,  block  thickness  and 
block  angle  of  attack.  We  also  believed  that  it  was 
necessary  to  characterize  the  fluid  dynamic  pressure 
by  calculating  both  the  average  pressure  magnitude 
and  center  of  pressure  on  the  block  bottom  sruface 
for  several  attack  angles,  from  horizontal  through 
the  angle  of  maximum  hydrostatic  rigjiting  mo¬ 
ment. 

Setup 

The  experiments  were  conducted  in  CRREL's 
warm  flume,  having  a  cross  section  of  0.91  by  0.91 
m,  a  total  length  of  7.32  m  and  a  variable  pumping 
capacity  up  to  0.3  m^/s.  A  hollow  Plexiglas  block 
"ice  floe,"  61.6  cm  square  by  7.62  cm  thick,  was 
constructed  with  91  pressure  taps  in  its  bottom 
surface.  The  taps  were  sized  to  accept  6.35-mm 
o.d.  polyethylene  manometer  tubing  and  were 
variably  spaced  as  shown  in  Figure  4.  A  Plexiglas 
cover  extended  over  the  top  front  half  to  keep  wa¬ 
ter  from  hitting  the  tubes  as  the  block  front  sub¬ 
merged  during  rotation.  The  block  was  fastened 
to  the  flume  by  hinged  supports  at  its  rear  and  an 
adjustable  threaded  rod  at  its  front  top  (Fig.  5). 
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Figure  4.  Tap  locations.  They  were 
variably  spaced  in  both  the  X  and  Y 
directions. 


Figure  5.  Block  attached  to  the  flume 
at  the  front  by  a  threaded  rod  and  at 
the  rear  by  hinged  supports.  These 
held  the  block  rigid  during  a  test  but 
still  allowed  the  attack  angle  to  be 
changed. 


This  allowed  both  vertical  and  rotational  (angle  of 
attack)  adjustments. 

Pressure  was  measured  by  routing  24  water- 
filled  manometer  tubes  via  a  fluid  switch  multi¬ 
plexer  to  a  liquid-to-gas  pressure  transducer  that 
had  a  ±6.35-cm  range.  In  this  manner  it  was  possible 


to  measure  the  pressure  of  24  separate  manometer 
taps  with  one  transducer.  A  stilling  tube  gave  us  a 
stable  measure  of  the  piezometric  head  in  the  flume, 
which  was  measured  at  the  beginning  and  end  of 
each  group  of  readings.  In  this  way  the  dynamic 
fluid  pressure  difference  from  static  was  determined 
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Table  1.  Test  matrix  showing  thickness-to-depth  ratios  (tlH),  velocities  (V)  and  angles  of  attack  used. 


t/H 


Velocity  V 
(m/s) 


Angle  of  attack 
(degrees) 


0.10 

0.45 

0.0 

3.3 

6,6 

10.0 

13.4 

16.8 

20.4 

0.13 

0.14* 

0.23 

0.39 

0.50 

0.0* 

3.3 

6.6 

10.0 

13.4 

16.8 

20.4 

0.20 

0.39** 

0.50 

0.60 

0.00** 

1.1 

3.3 

6,6 

10.0 

13.4 

16.8 

*  One  test  at  t/H  =  0.13,  V  =  0.14,  and  0°  angle  of  attack 
**  Four  tests  at  t/H  =  0.20,  V  =  0.39,  and  0°  angle  of  attack 


Figure  6.  Surface  plot  showing  a  typical  pressure  distribu¬ 
tion  across  the  full  width  of  our  block  during  a  symmetry 
check.  Only  the  first  three  rows  of  taps  were  used;  V  =  0.35  m/s, 
t/H  =  0.05  and  Q  =  0°. 

for  each  tap.  The  readings  and  their  standard  devia¬ 
tions  were  plotted  in  real-time  and  stored  on  a  per¬ 
sonal  computer. 

Test  procedure 

We  began  a  typical  test  by  setting  the  water  to  the 
desired  depth  and  horizontally  submerging  the 
block  to  the  required  displacement.  The  flume  ve¬ 
locity  was  adjusted  to  the  desired  value  and  two 
scans  of  25  pressures  (two  for  static  height  and  23 
taps)  were  made.  The  model  floe  was  rotated  to  the 
next  desired  angle  and  two  more  scans  made.  This 
procedure  was  repeated  until  aU  desired  block 
angles  were  examined  for  any  given  depth  and  ve¬ 
locity. 

We  did  the  experiments  in  two  phases.  The  first 
was  to  verify  the  symmetry  of  the  pressures  across 
the  fuU  width  of  the  block.  The  second  was  to  mea¬ 
sure  pressures  on  one  half  of  the  block,  assuming 


symmetry  of  the  pressure  distribution  for  the  other 
half. 


RESULTS 
Phase  one 

A  typical  pressure  distribution  obtained  during 
phase  one  of  our  experimental  program  is  shown  in 
Figure  6.  The  water  flow  was  in  the  positive  Y-direc- 
tion.  The  readings  were  obtained  at  a  zero  angle  of 
attack  and  a  velocity  of  0.35  m/s.  In  these  experi¬ 
ments,  a  negative  pressure  tended  to  rotate  the  block 
downward,  i.e.,  toward  the  negative  Z-direction;  a 
positive  pressure  pushed  against  the  block  in  the 
positive  Z-direction.  Examination  of  the  figure 
points  out  the  pressure  distribution's  high  degree  of 
symmetry  about  the  block  centerHne  at  X  =  30.5  cm. 
Having  proven  symmetry,  we  instrumented  one 
half  of  the  block  along  its  entire  length,  and  entered 
the  second  phase  of  the  experimental  program. 

Phase  two 

The  experiments  were  run  at  three  values  off/  H — 
0.10,  0.13  and  0.20 — and  were  obtained  from  our 
fixed  block  thickness  of  7.62  cm  and  three  different 
flow  depths.  Velocities  used  were  0.14,  0.23,  0.39, 
0.45, 0.50  and  0.60  m/s  (Table  1).  Block  attack  angles 
ranged  from  0°  to  approximately  20°.  Two  repUcates 
were  performed  at  each  condition,  except  as  noted 
in  Table  1,  and  the  data  presented  below  are  the  av¬ 
erage  of  these  two  replicates. 

Figure  7  fllustrates  a  typical  surface  pressure  plot 
obtained  during  phase  two.  Generally,  at  zero  angle 
of  attack,  the  pressures  were  lower  at  the  front  of  the 
block,  becoming  less  negative,  if  not  positive,  to¬ 
wards  the  rear  of  the  block.  TTie  rear  central  portion 
of  the  block  was  frequently  near  zero  pressure,  with 
perhaps  some  areas  of  negative  pressure  at  the  cor- 
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Figure  7.  Typical  pressure  distribution  surface  plot.  This  is 
created  by  reading  selected  taps  along  one  side  of  the  block  and 
generating  a  mirror  image  of  the  pressures  on  the  opposite 
side.  In  this  example  V  =  0.39  m/s,  t/H  =  0.13  and  0  =  0°. 

nets.  The  pressure  at  the  front  central  portion  of  the 
block  was  frequently  less  negative  than  at  the  front 
sides  or  the  block's  midsection.  As  the  angle  of  at¬ 
tack  increased,  the  magnitude  and  area  of  negative 
pressure  would  increase  and  move  towards  the  rear 
of  the  block  rmtU  eventually  the  entire  pressure  dis¬ 
tribution  was  negative.  Figure  8  shows  typical  pres¬ 
sure  distributions  from  our  tests:  those  from  test  1 
{t/H  =  0.10)  tended  to  be  more  uniform  between  the 
front  and  rear  sections  of  the  block  than  those  from 
tests  2  and  3  {t/H  =  0.13  and  0.20  respectively).  In 
test  3  the  pressure  difference  between  the  front  and 
rear  sections  of  the  block  is  greater  and  the  overall 
pressure  is  lower.  There  are  also  numerous  areas  of 
extremely  low  pressure,  probably  associated  with 
vortex  formation  and  flow  separation. 

To  characterize  the  pressure  distributions,  we  cal¬ 
culated  the  surface  area  over  which  the  pressure  acts 
and  the  average  for  regions  of  positive  and  negative 
pressure  (Fig.  9).  It  is  interesting  that  tests  1  and  3 
have  roughly  the  same  area  of  negative  pressure, 
with  a  much  greater  magnitude  for  test  3. 

Figure  10  shows  the  positive  and  negative  center 
of  pressure  (COP)  locations  (measured  from  the  ro¬ 
tation  axis)  and  the  resultant  moment  for  each  test. 
Here,  it  can  be  seen  that,  although  test  2  did  not  have 
quite  the  area  of  negative  pressure  that  test  1  did,  its 
negative  COP  is  farther  from  the  axis  of  rotation, 
while  its  positive  COP  is  about  the  same  as  test  I's. 
This  results  in  the  two  moments  being  very  similar 
in  the  initial  stages  of  rotation,  with  test  2's  being 
greater  at  the  higher  attack  angles  where  the  effect  of 
depth  is  probably  being  felt.  Test  3  has  a  much 
higher  moment  than  the  others,  but  its  negative 
COP  is  almost  identical  to  test  I's.  It  should  be  noted 


a.  Test  1:  t/H  =  0.10,  V  =  0.45  m/s. 

Figure  8.  Representative  pressure  distributions  from  three 
tests.  The  block's  trailing  edge  is  atY  =  0. 
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Pressi 


Figure  9.  Areas  of  positive  and  negative  pressure  on  the 
block's  lower  face,  made  nondimensional  by  the  total  block 
area.  Also  shovm  are  the  average  pressures,  made  non- 
dimensional  with  the  stagnation  pressure  of  the  flow. 

that  the  positive  COPs  in  tests  1  and  2  are  almost 
identical,  while  that  for  test  3  moves  away  from  the 
rotation  axis.  This  has  the  effect  of  hindering  the 
block's  rotation  if  we  assume  a  pinned  lower  comer. 

Since  we  assume  that  the  block  rotates  about  its 
downstream  lower  comer,  these  positive  and  nega¬ 
tive  pressures  cause  moments  about  that  point.  We 
therefore  use  our  moment  coefficient  a  as  an  indica¬ 
tion  of  the  total  moment  generated  by  the  pressures 
about  that  axis.  Figure  11  shows  a  plotted  against 
angle  of  attack.  There  is  an  obvious  trend  of  increas¬ 
ing  a  with  increasing  angle  of  attack.  This  trend  is 
exaggerated  in  shallow  water  (t/H  =  0.20)  but  is  still 
present  in  the  deeper  water.  There  also  appears  to  be 
a  weak  trend  of  decreasing  a  with  increasing  veloc¬ 
ity  at  any  fixed  t/ H  ratio.  Overall,  depth  seems  to 
have  a  larger  influence  upon  a  than  does  velocity, 
showing  a  trend  of  decreasing  a  with  increasing 
depth. 


Figure  10.  Center  of  negative  and  positive  pressure  di¬ 
vided  by  block  length  along  ivith  the  total  undertuming 
moment,  made  nondimensional  by  the  maximum  hydro¬ 
static  righting  moment  fCOP/L = 0  is  the  axis  cf  rotation). 

To  see  if  the  block  is  stable,  we  compare  the  mea¬ 
sured  moments  attributable  to  dynamic  fluid  pres¬ 
sure  with  the  block  hydrostatic  righting  moments  in 
Figure  12.  Both  the  measured  imdertuming  mo¬ 
ment  and  the  righting  moment  are  made  nondimen¬ 
sional  by  the  maximum  righting  moment  given  by 
eq  13.  Also,  the  data  are  presented  with  the  95%  con¬ 
fidence  interval,  which  quantifies  the  rmcertainty  in 
our  static  pressure  and  individual  tap  readings.  This 
uncertainty  will  be  discussed  below. 

It  can  be  seen  from  Figure  13  that  the  moments 
created  on  the  block  in  the  deeper  water  (f/  H  =  0.10 
and  t/H  =  0.13)  would  not  be  sufficient  to  overcome 
the  righting  moment.  At  zero  angle  of  attack,  the 
moment  would  start  the  block  rotating,  but  at  the 
next  measured  angle,  the  moment  is  below  the 
righting  moment  curve,  which  means  it  would  be 
stable.  Neglecting  rotational  inertia,  we  see  that  the 
block  is  stable  if  the  undertuming  moment  is  not  al- 
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ways  greater  than  the  righting  moment.  In  the  shal¬ 
low  water  {t/H  =  0.20),  however,  the  moments  gen¬ 
erated  at  every  velocity  of  our  experimental  pro¬ 
gram  would  be  sufficient  to  undertum  the  block. 

DISCUSSION 
Data  imcertainty 

Inspection  of  the  imcertainty  diown  in  Figure  12 
reveals  an  interesting  trend.  At  t/H  =  0.20  (shallow 
water),  the  imcertainty  is  generally  uniform  and 


1.5  I - ^ ^ ^  r 

b. 

t/H  =  0.13 

■  V  =  0.14  m/s 
1.0  —  •  0.23 

O  0.39 

O  0.50 


Angle  of  Attack 

Figure  11.  Nondimensioml  moment  co^cient  vs.  angle  of 
attack. 


-2  0  4°  8°  12°  16°  20° 

Angle  of  Attack 

a.  t/H=  0.13. 

Figure  12.  Nondimensional  undertuming  moment  vs. 
angle  of  attack. 
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small.  In  the  moderate  depth  water  oit/H  =  0.13,  a 
slightly  increasing  imcertainty  interval  is  evident 
with  increasing  velocity.  In  our  deepest  water  ex¬ 
periments  at  t/H  =  0.10,  Figure  12c  shows  that  the 
scatter  is  larger  and  takes  a  discontinuous  jump  at 
an  attack  angle  of  13.4°. 

We  believe  the  generally  larger  uncertainty  is 
caused  by  vibration  in  our  experimental  setup  at 
high  volumetric  flow  rates.  Our  block  was  attached 
to  a  frame  that  was  in  turn  attached  to  the  flume  side 


Angle  of  Attack 


walls.  At  the  deeper  water  condition  of  t/H  =  0.10 
and  V  =  0.45  m/s,  the  high-frequency  pump  vibra¬ 
tion  was  transmitted  by  the  flume  piping  to  the 
flume  walls  and  thus  to  our  block.  The  vibration 
was  periodic  and  manifested  itself  as  an  increase  in 
the  calculated  standard  deviation  of  the  differential 
pressures.  We  believe  that  the  accuracy  of  the  aver¬ 
age  values  was  not  affected  because  we  were  able  to 
analyze  the  readings  in  real-time  during  each  ex¬ 
periment  to  ensure  that  the  sampling  rate  was  ad¬ 
equate  for  the  vibration  frequency.  Our  replicate  val¬ 
ues  used  to  obtain  the  average  shown  in  the  figures 
are  close  to  each  other.  In  the  worst-case  deep  water 
tests,  the  average  difference  between  the  replicates 
and  the  averaged  values  was  39%,  ranging 

down  to  less  than  10%  for  the  shallow  water  tests, 
giving  us  confidence  in  the  procedure. 

Static  analysis  of 
block  undertuming 

Using  the  above  results,  we  look  at  block 
undertuming  through  a  static  moment  analysis,  as¬ 
suming  a  clockwise  rotation  as  positive,  and  know¬ 
ing  that  for  a  rotating  block 

IM>0  (15) 

where  M  is  a  moment  acting  on  the  block.  We  ex¬ 
press  these  moments  as  a  function  of  block  angle  of 
attack  0  by 

Mf(0)-Mr(e)>O  (16) 

with  Mf  being  the  imdertuming  moment  created  by 
the  fluid  flow  and  Mj.  the  hydrostatic  righting  mo¬ 
ment.  By  neglecting  the  inertial  component,  we 
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b.  t/H  =  0.20. 


c.  t/H  =  0.10. 


Figure  12  (cant'd). 
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Figure  13.  Nondimensional  under- 
turning  moments  vs.  angle  of  attack 
for  different  t/H  ratios  at  similar  ve¬ 
locities. 


Figure  14.  Froude  criterion  for  block 
stability. 


have  a  block  instability  criterion  for  all  angles  that 
will  under-predict  when  a  block  becomes  unstable, 
that  is 


and 


M,(e»)= 


(20) 


.^r(9)J 


>1 . 


(17) 


By  inspection  of  the  data,  we  can  see  that  the 
tmdertuming  moment  is  Unear  through  the  angle 
where  the  hydrostatic  righting  moment  is  a  maxi¬ 
mum.  This  allows  us  to  restrict  eq  17  to  0*,  the  angle 
of  maximum  righting  moment,  or 


'Mf(9*)' 


>1  . 


(18) 


With  substitution  of  eq  2, 5, 6,  and  13  into  eq  18,  we 
arrive  at 


Mf(9*): 


ic(0*)p^y2 


(71)[COPy(0=^)]  (19) 


where  pj  is  ice  density  and  g  is  the  acceleration  of 
gravity.  The  data  indicate  that  for  a  particular  t/H, 
both  C  and  COPy  are  functions  of  the  angle  of  attack 
9.  Taking  this  into  account  and  substituting  eq  19 
and  20  into  18,  we  arrive  at  the  expression 


r 


A 


y 


f 

1- 

Pi 

iPw  J 

y 

V 

C(9*)(COPY(9^)j^^  ^ 


(21) 
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which  in  the  form  of  a  densimetric  Froude  number 
(Fj)  is  as  follows 


fja(e’)]>i. 


(22) 

And  by  defining 

f/= - (23) 

(a(e*))2 

we  know  that  for  block  stability  Fj.  must  be  less  than 
fj.*  at  any  given  t/H  value.  In  Figure  14,  is  plot¬ 
ted  against  at  t/H  =  0.13  and  t/H  =  0.20  along 
with  a  line  fit  to  the  points  by  linear  regression.  Our 
one  data  point  at  t/H  =  0.10  is  included  but  was  not 
used  to  find  the  regressed  line.  The  45°  line  depicts 
the  critical  Froude  number  for  stability.  The  effect  of 
the  depth  change  onf,*  between  t/H  =  0.20  and  t/H 
=  0.13  is  very  evident,  while  increasing  the  depth 
from  t/H  =  0.13  to  t/H  -  0.10  would  appear  to  have 
no  effect.  At  t/H  -  0.20,  F^*  is  less  sensitive  to  veloc¬ 
ity  than  at  t/H  =  0.13.  We  would  expect  that  more 


data  at  t/H  =  0.10  would  show  little,  if  any,  change 
from  t/H  =  0.13. 

Dynamic  analysis  of 
block  undertuming 

Figure  15  shows  the  moment  coefficient  as  a 
function  of  block  angle  of  attack  0  for  three  thick- 
ness-to-depth  ratios  and  several  velocities.  By  ne¬ 
glecting  the  velocity  dependence,  the  moment  coef¬ 
ficient  becomes  a  function  of  angle  of  attack  and 
thickness-to-depth  ratio  (Fig.  16).  A  cubic  polyno¬ 
mial  fit  of  the  moment  coefficient  was  made  to  rep¬ 
resent  the  data  as  a  continuous  function  of  angle  of 
attack  (in  degrees).  The  curve  fit  takes  the  form 

—  =  Ao+Ai0  +  A2e2+A303  (24) 

with  the  coefficients  shown  in  Table  2. 


Angle  of  Attack  (degrees) 


Figure  15.  Moment  coefficient  for  two 
thickness-to-depth  ratios  and  several 
velocities. 


0  5°  10°  15°  20°  25° 

Angle  of  Attack 
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Figure  16.  Cubic  polynomial  curve  fit 
of  moment  coefficients  to  block  angle  of 
attack. 


Table  2.  Cubic  polynomial  curve  fit  of  moment  coefficient  with 
angle  of  attack  in  degrees. 


t/H 

Aq 

Ai 

A2 

A3 

R2 

0.10 

2.33x10-1  , 

-1.63x10-2 

3.75x10-3 

-1.24x10-1 

0.922 

0.13 

1.96x10-1 

9.56x10-3 

4.73x10-^ 

8.48x10-2 

0.994 

0.20 

2.04x10-1 

8.67x10-2 

-4.49x10-3 

1.90x10-1 

0.998 

Moment  analysis 

The  equation  for  the  motion  of  the  block  is 

lM  =  l'd  (25) 

where  the  moment  of  inertia  I  is  given  by 

1  =  |pi  LA  (L2+f2)  (26) 

subject  to  the  initial  conditions 

00  =  0 

0  =  0 

q_Mu(0o)-M,(0o)  ^27) 

In  the  expression  for  the  moment  of  inertia,  we 
have  made  tiae  assumption  that  fluid  inertia  may  be 
neglected.  Our  pressure  measurements  quantified 
the  acceleration  of  the  fluid  around  the  block;  there¬ 
fore,  we  are  neglecting  only  the  fluid  acceleration 
associated  with  block  angular  velocity  and  accelera¬ 
tion.  This  is  reasonable  because  the  upstream  flow 
velocity  was  estimated  to  be  approximately  five 
times  greater  than  the  maximum  velocity  of  the 
block  leading  edge  as  it  rotates,  which  means  that  the 
fluid  acceleration  through  time  is  relatively  small 
compared  to  the  acceleration  in  space.  An  estimate 
of  the  fluid  added  mass  would  require  a  knowledge 
of  the  volume  of  recirculating  flow,  which  changes 
with  both  time  and  angle  of  attack.  Such  an  estimate 
would  be  challenging  while  giving  only  marginal 
improvement  in  the  solution  of  eq  25. 

It  is  useful  to  define  the  external  work  done  on  the 
block  by  the  fluid.  If  a  mechanical  energy  balance  is 
done  on  an  ice  floe  that  is  initially  at  rest,  we  find  that 
the  cumulative  work  done  on  the  block  is  equal  to 
the  rotational  energy  of  the  block.  Formally,  this  is 
expressed  as 

i702=J(Mu-Mr)a0.  (28) 

2  0 


Evaluation  of  eq  28  leads  to  the  five  char¬ 
acteristic  energy  curves  shown  in  Figure  17. 
Each  kinetic  energy  curve  is  shown  with  its 
corresponding  righting  and  undertuming 
moments.  Figure  17a  shows  that  the  under¬ 
tuming  moment  is  greater  than  the  righting 
moment  for  all  angles  considered.  This  con¬ 
dition  leads  to  unequivocal  block  instabil¬ 
ity,  as  is  evidenced  by  the  continuously  in¬ 
creasing  kinetic  energy  curve. 

Upon  lowering  the  velocity,  the  imdertuming  mo¬ 
ment  becomes  tangent  to  the  righting  moment  at 
exactly  one  point,  as  is  shown  in  Figure  17b.  At  the 
point  of  tangency,  the  kinetic  energy  curve  achieves 
a  point  of  inflection  where  the  angular  acceleration 
equals  0.  In  the  absence  of  rotational  inertia,  the 
block  would  remain  at  this  angle  of  attack,  since  the 
undertuming  moment  equals  the  righting  moment. 
A  somewhat  lower  velocity  leads  to  a  righting  mo¬ 
ment  in  excess  of  the  undertuming  moment,  which 
means  that  the  block  would  be  stable.  A  higher  ve¬ 
locity  would  lead  to  a  larger  undertuming  moment, 
which  means  that  the  block  would  be  unstable.  There¬ 
fore,  this  represents  the  point  of  static  metastability. 

Figure  17c  shows  the  effect  of  a  further  reduction 
in  approach  velocity.  In  this  case,  the  block  achieves 
a  local  maximum  in  kinetic  energy  at  the  first  in¬ 
tersection  of  undertuming  moment  and  righting 
moment.  Then,  the  righting  moment  is  greater  than 
the  undertuming  moment,  leading  to  a  slowing  of 
the  block,  as  is  evidenced  by  a  reduction  in  block 
kinetic  energy.  Finally,  at  the  second  intersection  of 
the  righting  and  undertuming  moments,  the  kinetic 
energy  achieves  a  local  minimum.  Beyond  this  point, 
the  block  is  unstable.  Clearly,  it  is  the  block  rotational 
inertia  that  carries  the  block  through  these  angles 
where  the  righting  moment  exceeds  the  undertuming 
moment. 

If  we  reduce  the  velocity  still  further  (Fig.  17d), 
the  block  again  achieves  a  local  maximum  in  kinetic 
energy,  but  this  time  the  local  minimum  occurs  at  a 
kinetic  energy  of  0.  This  means  that  the  block  has 
been  brought  to  rest  at  a  positive  angle  of  attack 
where  the  undertuming  moment  equals  the  right¬ 
ing  moment.  This  point  represents  a  metastable  po¬ 
sition,  since  a  perturbation  in  velocity  will  either 
render  the  block  stable  or  xmstable.  Therefore,  this 
velocity  and  angle  of  attack  may  be  considered  criti¬ 
cal  to  dynamic  block  stability. 

Finally,  we  reduce  the  velocity  to  achieve  the 
stable  condition  shown  in  Figure  17e.  In  this  case, 
the  block  picks  up  angular  momentum  and  achieves 
a  local  maximum  in  kinetic  energy  as  before.  Flow- 
ever,  the  kinetic  energy  of  the  system  never  achieves 
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Figure  17.  Characteristic  energy  curves. 


a  local  minimum.  Instead,  it  becomes  equal  to  zero 
at  an  angle  where  the  righting  moment  exceeds  the 
undertuming  moment.  This  condition  can  cause  the 
block  to  oscillate  stably  in  this  range  of  angular  dis¬ 
placement. 

To  evaluate  the  effect  of  rotational  inertia  on 
block  stability,  we  define  the  densimetric  Froude 
number  at  metastability  f  /  as 


Comparison  with 
previously  published  results 

Other  investigators  have  given  block  stability  cri¬ 
teria  that  may  be  compared  with  the  current  analy¬ 
sis.  Pariset  and  Hausser  (1961)  gave  the  following 
stability  criterion  in  the  form  of  a  Froude  number 
with  depth  as  the  characteristic  length  scale.  Their 
equation  is 


fr 


V 


(29) 


1-- 


Pw 


with  V*  taking  on  a  value  associated  with  Figure  17b  for 
static  stability  criteria  and  Figure  17e  for  dynamic 
stability  criteria.  An  iterative  procedure  was  imple¬ 
mented  to  determine  the  velocity  that  gave  the  proper 
conditions  for  both  static  and  dynamic  stability. 


V 


(30) 


where  H  -  upstream  depth  of  flow 
V  =  approach  velocity 
K  =  empirical  parameter. 
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The  parameter  k  is  not  used  in  a  later  paper  by 
Pariset  et  al.  (1966). 


Kinetic  Energy  (J)  Moment  (N-m) 


Figure  17  (cant’d).  Characteristic  energ]/ curves. 


Ashton's  stability  criterion  is  in  the  form  of  a  den- 
simetric  Fronde  number 


ting  constant  equal  to  0.27.  For  deep  water  condi¬ 
tions,  Daly  used 


y 


r  \ 

f'l 

i_2l 

1  Pw  7 

(31) 


Daly  (1984)  introduced  a  shallow-water  stability 
criterion  of  the  form 


where  t'  is  block  freeboard  and  P  is  an  empirical  fit- 


(33) 


where  Sf  is  the  Strouhal  number  at  block  instability, 
and  may  be  determined  from  the  data. 

In  Coutermarsh  (1986),  a  value  of  2.08  was  as¬ 
sumed  for  a  wide  variety  of  block  geometries,  while 
in  the  current  work  a  value  of  1.78  is  used.  Figure 
18  shows  these  stability  criteria  and  the  critical 
Froude  numbers  determined  from  the  current 
analysis.  As  can  be  seen,  both  Pariset  and  Hausser 's 
and  Ashton's  stability  criteria  under-predict  the 
critical  Froude  number  for  the  block  geometry 
under  consideration. 
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Kinetic  Energy  (J)  Moment  (N-m) 


e.  Stable  block  at  t/H  =  0.13  and  V  =  0.48  m/s. 
Figure  17  (cant'd). 


Daly's  criterion  does  reasonably  well  in  predict¬ 
ing  both  the  deeper-water  and  the  shaUow-water 
cases.  However,  if  Daly's  deep-water  criterion  is 
used  with  a  Strouhal  number  of  2.08,  it  also  imder- 
predicts  the  critical  Froude  numbers.  Because  in  the 
current  work  we  neglect  the  effects  of  hydrody¬ 
namic  fluid  pressure  on  the  top  of  the  block  and  on 
the  front  face,  the  block  is  expected  to  undertum  at 
slightly  lower  Froude  numbers  than  those  calcu¬ 
lated.  Indeed,  fluid  pressure  on  the  top  of  the  block 
would  reduce  the  difference  between  the  calculated 
critical  Froude  number  and  Daly's  predictions 
based  on  a  Strouhal  number  of  2.08.  In  any  case,  the 
Strouhal  number  is  a  function  of  block  geometry, 
and  a  variation  between  2.08  and  1.78  is  possible, 
given  the  differences  in  block  geometry.  Tlie  most 
salient  feature  of  tire  deep  water  data  is  an  apparent 
leveling  off  of  the  critical  Froude  number  with  de¬ 
creasing  thickness-to-depth  ratio.  This  has  an  intui¬ 
tive  appeal  since  we  expect  that  the  moment  coeffi¬ 
cient  should  become  completely  insensitive  to  arbi¬ 
trary  increases  in  fluid  depth  for  very  deep  water 
conditions.  The  current  work,  coupled  with  Daly's 
deep-water  criterion,  lends  support  to  this  argu¬ 
ment. 

It  is  difficult  to  know  exactly  how  weU  Daly's  cri¬ 
terion  represents  the  shaUow-water  data  since  we 
have  only  one  such  data  point.  In  spite  of  this  draw¬ 
back,  Daly's  criterion  does  appear  to  represent  the 
expected  trend  reasonably  well.  More  data  are  re¬ 
quired  before  a  definite  conclusion  may  be  drawn 
about  his  sKallow-water  criterion. 

An  important  feature  of  Figure  18  is  the  offset  be¬ 
tween  the  static  and  dynamic  Froude  numbers  re- 


Figure  18.  Comparison  of  calculated  stability  points  mith  existing  stability 
criteria. 
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these  angles  would  hold  for  other  block  ge¬ 
ometries,  it  appears  that  the  upper  limit  is 
not  far  from  the  angle  at  which  the  maxi¬ 
mum  righting  moment  occurs. 

CONCLUSIONS 

The  hydrodynamic  fluid  pressures  be¬ 
neath  an  idealized,  parallelepiped  model  ice 
block  were  measured.  The  ice  block  was 
fixed  in  space,  thereby  eliminating  angular 
acceleration  during  measurement  of  the 
fluid  pressures.  The  pressure  distribution 
was  shown  to  be  symmetrical  about  the 
Thickness  to  Depth  (t/H)  centerline  of  the  block  and  included  regions 

turning  moment  for  the  entire  pressure  dis- 
quired  to  destabilize  the  block.  Evidently,  for  these  tribution  was  compared  with  the  hydrostatic  right- 

thickness-to-depth  rahos  and  block  geometry,  fluid  ing  moment,  and  the  condition  of  static  stability  was 

inertia  is  responsible  for  a  5  to  10%  decrease  in  the  described.  A  flow  area  reduction,  which  includes 

Froude  number  to  destabilize  the  block.  Further-  the  thickness-to-depth  ratio,  the  length-to-depth  ra- 

more,  we  expect  a  further  5  to  10%  decrease  in  criti-  tio  and  the  angle  of  attack,  is  presented.  A  moment 

cal  Froude  number  by  includmg  the  hydrodynamic  coefficient  was  defined  that  relates  the  magnitude 

fluid  forces  on  the  top  and  front  face  of  the  block,  and  point  of  application  of  the  dynamic  fluid  pres- 

which  have  been  ignored  in  the  current  work.  sures.  It  is  evaluated  in  relation  to  flow  depth  and 

Figure  19  shows  the  angle  at  which  metastability  flow  velocity  for  a  fixed  block  geometry.  The  mo- 

occurs.  This  angle  is  compared  to  other  often-cited  ment  coefficient  appears  to  have  a  linear  relation- 

angles  such  as  the  no-spill  angle  and  the  angle  of  ship  with  the  reduction  ratio  over  the  range  of  pa- 

maximum  righting  moment.  The  angle  of  static  rameters  investigated.  This  linear  relationship  sug- 

metastability  is  lower  than  the  angle  of  dynamic  gests  that  block  length  is  as  important  as  thickness 

metastability  because  of  how  each  is  defined.  Static  for  characterizing  block  stability  in  deep  water  con- 

metastability  occurs  at  the  point  of  tangency  of  the  ditions,  especially  for  high  1/ 1  aspect  ratio  blocks, 

righting  and  imdertuming  moments.  E^amic  The  imdertuming  phenomenon  is  related  to  the 

metastability  occurs  at  the  second  intersection  of  described  hydrostatic  righting  moment  of  the  block, 

righting  and  imdertuming  moments  that  gives  zero  which  depends  on  block  geometry  and  density  and 
block  kinetic  energy.  These  conditions  are  shown  fluid  density.  By  combining  the  two  phenomena,  we 
graphically  in  Figures  17b  and  d.  arrive  at  a  block  underturrring  description  in  terms 

It  is  interesting  to  note  that  at  deep  water  condi-  of  moments.  We  have  shown  that  this  block  un- 

tions  the  angle  of  maximum  righting  moment  may  dertuming  can  be  described  by  a  relatively  simple 

be  very  close  to  the  angle  of  dynamic  metastability,  equation  that  relates  a  Froude  number  to  block  in¬ 
while  at  very  shallow  water  conditions  the  classic  stability.  Furthermore,  if  the  block  were  not  con- 

no-spill  angle  may  be  very  close  to  the  metastability  strained  at  its  lower  downstream  comer,  the  couple 

angle.  More  data  are  required  before  these  relation-  developed  by  the  fluid  forces  would  tend  to  rotate  it 

ships  may  be  established.  with  less  effort  by  uplifting  the  rear  of  the  block. 

The  angle  at  which  metastabUity  occurs  may  This,  along  with  our  not  including  any  rotational  in- 

have  little  practical  importance  except  to  guide  fur-  ertia  effects  in  the  above,  means  we  have  a  moment 

ther  research  that  focuses  on  measuring  the  un-  analysis  that  will  somewhat  under-predict  when  a 

dertuming-righting  moments  at  various  angles  of  block  will  undertum. 

attack.  For  example,  it  appears  that  for  this  block-  Our  dynamic  analysis  of  stability  shows  that  ice 

length-to-thickness  aspect  ratio  and  density,  angles  floe  angular  momentum  may  lead  to  rotational  in- 

greater  than  10°  need  not  be  considered  in  an  experi-  stability  even  when  the  righting  moment  is  greater 

mental  program.  Although  it  is  unknown  how  well  than  the  undertuming  moment  at  a  given  angle  of 
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attack.  In  this  analysis,  the  net  work  on  the  block 
was  equated  to  the  rotational  kinetic  energy  of  the 
block.  This  yielded  estimates  of  stability  for  the 
model  ice  floe  for  three  thickness-to-depth  ratios, 
which  were  shown  to  compare  favorably  to  an  exist¬ 
ing  stability  criteria.  The  angular  momentum  causes 
block  instability  at  Froude  numbers  5  to  10%  below 
those  associated  with  static  stability.  The  calculated 
Froude  numbers  are  assumed  to  over-predict  the  ac¬ 
tual  critical  Froude  numbers  since  the  hydrody¬ 
namic  pressure  distribution  on  the  top  and  front  of 
the  model  ice  floe  was  not  measured. 

The  measured  areas  of  positive  pressure  at  the 
back  of  the  block  suggest  that  it  is  possible  for  a 
block  to  be  pushed  out  of  the  water  vertically  at  its 
downstream  edge.  Researchers*  we  have  spoken 
with  have  confirmed  seeing  this  phenomenon. 

We  were  constrained  in  our  ability  to  present  a 
complete  set  of  moment  coefficients  by  our  lim¬ 
ited  experiments.  Future  work  in  this  area  would 
do  well  to  investigate  different  block  thickness-to- 
depth  ratios  and  block  lengths.  The  three-dimen¬ 
sional  pressure  distributions  also  suggest  that  the 
effect  of  block  width  should  be  investigated.  Con¬ 
sideration  of  blocks  with  rounded  leading  edges 
might  also  prove  interesting  to  see  the  variation  in 
flow  attachment.  Better  flow  visualization  may 
also  shed  light  on  the  phenomenon  of  positive 
pressures  beneath  the  ice  block. 

Deep  water  data  (t/H  >  0.10)  at  velocities  over 
0.45  m/s  would  be  interesting  and  should  be  ob¬ 
tained.  We  were  restricted  in  both  available  depth 
and  achievable  velocities  in  our  facilities.  We,  there¬ 
fore,  could  not  fully  characterize  the  deeper  flows. 

As  mentioned  before,  our  experiments  did  not 
include  any  effect  attributable  to  the  dynamic  fluid 
impact  against  the  front  and  top  of  the  block.  Experi¬ 
ments  that  measure  force  on  the  block  at  a  known 
distance  from  the  assumed  rotation  point  would 
contain  this  information,  although  they  would  not 
reveal  the  nature  of  the  forces. 


*  Personal  communication  with  Dr.  G.  Ashton  and  Dr.  J. 
Lever,  both  of  CRREL,  1991. 


Because  one  of  the  ultimate  goals  of  this  work  is 
to  estimate  the  capture  efficiency  of  an  ice  cover,  fu¬ 
ture  efforts  should  also  be  focused  on  characterizing 
the  typical  ice  floe  geometries  formd  in  the  field  and 
then  quantifying  the  righting  and  undertuming  mo¬ 
ments  associated  with  these  geometries. 

Additionally,  the  initial  impact  of  an  ice  floe  with 
an  ice  cover  may  be  destabflizing  and  should  be  in¬ 
vestigated  as  a  separate  issue.  With  all  of  the  desired 
data,  a  family  of  curves  could  be  developed  to  more 
fully  document  the  capture  efficiencies  for  a  distri¬ 
bution  of  geometries  and  range  of  flow  conditions. 
Such  calculations  would  be  useful  for  designers  and 
operators  of  ice  control  structures,  as  well  as  for  nu¬ 
merically  modeling  the  behavior  of  floating  ice  floes. 
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